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Abstract. We examine the isentropic QCD equation of state within a quasi- 
particle model being adjusted to first principle QCD calculations of two quark 
flavours. In particular, we compare with Taylor expansion coefficients of en- 
ergy and entropy densities and with the isentropic trajectories describing the 
hydrodynamical expansion of a heavy-ion collision fireball. 
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1. Introduction 

With growing evidence, ideal hydrodynamics appears to be successfully describing 
the expansion stage of strongly interacting matter in relativistic heavy-ion collisions 
[QJ El El EJ El- The heart of hydrodynamics is the equation of state (EoS) 
which is needed as important interrelation among the state variables of strongly 
interacting matter in order to solve the hydrodynamic equations of motion. The 
EoS can, for instance, relate pressure, energy density and baryon density in the form 
p = p(e, Ub)- The influence of different model EoS on observables was analyzed, 
e. g. in [El El 0- Recent progress in QCD calculations (performed numerically on a 
discretized space-time grid, dubbed lattice), though, allows the calculation of these 
EoS quantities from first principles. It is the aim of the present paper to compare 
our quasi-particle model with some available lattice QCD results. 

Our paper is organized as follows. In section El we review our quasi-particle 
model. In section El we compare the model with recent two flavour lattice QCD 
results for pressure, energy density and entropy density. In particular, we focus on 
the isentropic trajectories and the EoS along those paths of constant entropy per 
baryon. The results are summarized in section El 
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2. Quasi-Particle Model 

The QCD EoS was often formulated in terms of quasi-particles with effectively 
modified properties due to the strong interaction (cf. [0 E] and references therein). 
More recently developed quasi-particle models are proposed in [^3 HH 1121 1131 1141 
1151 I16j . In our quasi-particle model (QPM) we employ as thermodynamic 

potential the pressure p in thermal equilibrium. Concentrating on the case of Nf = 2 
light quark flavours with one chemical potential (j, q , p as a function of temperature 
T and [i q reads 

V{T, fig) = ^ P a - B(T, fig) (1) 

a=q,g 

with partial pressures of quarks (q) and transverse gluons (g) 
Pa = d a /(6ir 2 ) J™ dfcfc 4 (/+ + /-) /w„. Here, d q = 2N f N c , d g = N 2 - 1, N c = 3, 
and = (exp([w a =p Ha]/T) + if? a ) -1 with S q = 1 for fermions and S g = —1 for 
bosons. Note \i g = 0. 

The quasi-particles propagate predominantly on-shcll with a dispersion rela- 
tion approximated by the asymptotic mass shell expression near the light cone, 
ijj a = \/k 2 + rria, where m a = m^. a + H a [ 118] with self-energies ±l a and mo- g = 0. 
For n a , we employ the asymptotic expressions of the gauge independent hard- 
thermal loop / hard-density loop self-energies [El- The mean field interaction 
term B(T, fi q ) in Q is determined by thermodynamic self-consistency and station- 
arity of the thermodynamic potential under functional variation with respect to the 
self-energies, 5p/8Tl a = [EJ. 

Replacing the running coupling g 2 entering II a by an effective coupling G 2 de- 
pending on T and fi q , non-perturbative effects within the strongly interacting system 
are accommodated in our model. A convenient parametrization of G 2 (T, \i q = 0) 
(cf. [Ell EU) is 

f Ga-joapKCr))' T * T - 
G 2 (T,^ = 0)= (2) 

I G|_ loop (e(T c )) + 6(1 - T/T c ), T < T c , 

where G2_ loop is the relevant part of the 2-loop running coupling. Here, £,(T) — 
\(T—T S )/T C contains a scale parameter A and an infrared regulator T s . The effective 
coupling G 2 for arbitrary T and fi q can be found by solving a quasi-linear partial 
differential equation which follows from Maxwell's relation (cf. [ 1171 123] for details 
of the model). In the next section, we want to test our phenomenological QPM by 
comparing with recent two-flavour (Nf = 2) lattice QCD results [1231 1251 12^|. 

3. Equation of State for N f = 2 

For small fj, q , the pressure can be expanded into a Taylor series in powers of (fx q /T), 

oo 

KT, M9 )=T 4 £^(T)(^) n , (3) 

n=0 
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Fig. 1. Comparison of our QPM with lattice QCD results (symbols) for cq(T) 
as function of T/T c for Nf — 2. Raw lattice QCD data [ 123 are continuum 
extrapolated by an extrapolation factor d = 1.1 for T > T c as advocated in [ 1251 127) 
due to finite size and cut-off effects. QPM parameters: A = 4.4, T s = 0.67T C , 
b = 344.4, B(T C ) = 0.31T C 4 , setting T c = 175 MeV as in [03]. The horizontal 
line depicts the corresponding Stefan-Boltzmann value highlighting the effects of 
strong interaction near T c . 

which has recently been studied in lattice QCD [ 1251 126| . The Taylor series was cal- 
culated up to including order (fi q /T) 6 . c„(T), vanishing for odd n, follow straight- 
forwardly from Q through differentiation, 



1 d n ( P /T 4 



n\ d( N /Ty 



(4) 



In Fig. ^ we compare QPM with lattice QCD results for cq(T) = p(T,fi q = 
0)/T 4 . In analogy to the lattice simulations, we set m,Q- q (T) = OAT. Furthermore, 
in [|21] an impressively good agreement between the QPM results of C2^fi{T) and 
the lattice QCD data was shown. 

From @, other thermodynamic quantities such as net baryon density tib = 
dp/d^,B, entropy density s and energy density e follow as 



s(T, MB )=T 3 ^. 



n=0 



l( T)(p n 



-3T/ 



eCZ> B )=T 4 ]T. 



n=0 



.(T,(i)" 



(5) 



where e„(T) = 3c n {T) + c' n (T) and s„(T) = (4 - n)c n {T) + c' n {T). Here, fi B = 3^ 9 
denotes the baryo-chemical potential and c' n (T) — Tdc n (T)/dT. Since s n (T) and 
e„(T) contain both, c n (T) and c^(T), they serve for a more sensitive test of the 
model than considering c„(T) alone. Estimating c' n (T) through a fine but finite 
difference approximation of c n (T), we compare QPM with lattice QCD results [ 
124) for S2,4 and e2,4 in Fig. [21 and find a fairly good agreement. The pronounced 
structures in the vicinity of the transition temperature follow from the change in 
the curvature of G 2 (T, \x q = 0) at T = T c . 
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T/T c 



Fig. 2. Comparison of our QPM with lattice QCD results for e„(T) (squares) 
and s n (T) (circles) as function of T/T c for Nf — 2; left (right) panel for n = 2 (4). 
QPM parameters: A = 12, T s = 0.87T C , b = 426.1 with T c = 175 MeV as adjusted 
to c 2 (T) from [ 26 (cf. [E2J). We choose AT = 5.25 MeV for suitably estimating 
c' n (T) by finite difference approximation of c n (T). 

Assuming local entropy and baryon number conservation during the hydro- 
dynamical expansion of the fireball created in heavy-ion collisions, the strongly 
interacting system evolves isentropically. The evolutionary paths of individual fluid 
elements can be displayed by trajectories s/ns — const in the T - fis plane. We 
calculate ns and s from (J5J up to 0(([Ib/T) 6 ) for the isentropic trajectories s/ub = 
300, 45 and compare with lattice QCD results in Fig. [21 Finding a fairly good 
agreement, the pattern of the evolutionary paths is mainly influenced by sq(T). For 
instance, in the vicinity of T c a larger value of so (T) of about 30% translates into a 
28% increase in [Ib for the same trajectory. At large T, where cq^{T) are essentially 
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Fig. 3. Isentropic trajectories: Comparison of Nf = 2 lattice QCD results [ I24| 
for s/tib — 300 (triangles) and 45 (circles) with the corresponding QPM results. 
The physical scale is T c = 175 MeV. 
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flat, the relation ^P- = holds for small /j,b, i- e. lines of constant specific 

entropy are given by lines of constant fis/T. 

Along the isentropic trajectories, we evaluate the EoS p(e). As depicted in 
Fig. 0] (left panel), the EoS, reproducing the lattice QCD results impressively well, 
is found to be almost independent of the considered specific entropy. Accordingly, 
the speed of sound = dp/de as exhibited in the right panel of Fig.0]is also rather 
independent of s/ub- 

4. Conclusions 

We presented a comparison of our QPM with recent two-flavour lattice QCD results 
of the isentropic equation of state at finite baryo-chemical potential. In particular, 
we focused on the Taylor expansion coefficients of energy density and entropy den- 
sity, reproducing the pronounced structures in the vicinity of the (pseudo-) critical 
temperature fairly well. In addition, the isentropic trajectories in the T - ^b plane 
were compared. The EoS along those paths of constant entropy per baryon was 
found to be rather independent of the particular value of specific entropy. Having 
tested the successful applicability of our model in the finite /is - region for Nf = 2, 
one can proceed to the physically interesting case of Nf = 2 + 1 quark flavours and 
extend the EoS p(T) towards finite baryon densities. Such an EoS can be applied 
to the hydrodynamical stage of heavy-ion collisions, e. g. by studying transverse 
momentum spectra and differential elliptic flow of various hadron species. This will 
be reported elsewhere. 




Fig. 4. Left panel: comparison of lattice QCD results [E3| of p as function of e 
for Nf = 2 along s/ns = 300 (triangles) and 45 (circles) with corresponding QPM 
results depicted by solid and dashed lines, respectively. Right panel: according 
speed of sound v 2 s . 
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